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Abstract. We consider a class of steady-state 
self-gravitating accretion disks for which efficient cooling 
mechanisms are assumed to operate so that the disk is self- 
regulated at a condition of approximate marginal Jeans 
stability. In an earlier paper, this scenario had been shown 
to lead naturally, in the absence of a central point mass, 
to a self-similar solution characterized by a flat rotation 
curve. In this article we investigate the entire parame- 
ter space available for such self-regulated accretion disks 
and provide two non-trivial extensions of the model. The 
first extension is that of a bimodal disk, obtained by par- 
tially relaxing the self-regulation constraint, so that full 
matching with an inner "standard" Keplerian accretion 
disk takes place. The second extension is the construc- 
tion of self-regulated accretion disks embedded in a diffuse 
spherical "halo". The analysis is further strengthened by 
a careful discussion of the vertical structure of the disk, in 
such a way that the transition from self-gravity dominated 
to non-gravitating disks is covered uniformly. 
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1. Introduction 

For simplicity, traditional models of accretion disks ignore 
the self-gravity associated with the accreting material and 
consider the case where the disk is approximately Keple- 
rian (e.g., see Pringle |1981| ). On the other hand, several 
studies have been carried out in order to investigate the 
effects related to the disk self-gravity. Studies of this type 
have addressed three different levels. The first stage is that 
where the disk self-gravity is incorporated in the study 
of the vertical structure of the disk. Here we may recall , 
among others, the pioneering work by Paczy hski (1978 ) 
and the very recent analysis by Bardou et al. ( 1998 ). The 
second important level is that where one considers the 
effects of the disk self-gravity on waves and on transport 
processes in the disk. Here we may recall a number of inter- 
esting analyses, in particular the studies by Lin & Pringle 
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cosity prescription proposed by Lin & Pringle (1990) (see 
also Andalib et al. 1997 and references therein). Recent hy- 
drodynamical simulations (Laughlin & Bodenheimer 1994, 
Laughlin & Rozyczka 1996| , Laughlin et al. 1997) have fo- 
cused on one key aspect of the modeling of accretion disks, 
i.e. they have tried to assess whether the standard con- 



stant a viscosity prescription (Shakura & Sunyaev 1973) 



is justified in systems where non-axisymmetric instabili- 
ties driven by the disk self-gravity play a major role. The 
third level where the self-gravity of the disk can be con- 
sidered is that where, because of the gravitational field 
contributed by the disk matter, the rotation itself of the 
accreting material is affected and is no longer kept to be 
Keplerian. In this direction ve ry few attempts have been 
made (see, e.g., Bodo & Curir 1992). 

The need for investigations where the self-gravity of 
the disk is properly taken into account is further stimu- 
lated by recent accurate observations that point to signifi- 
cant deviations from Keplerian rotation in objects that are 
naturally interpreted as accretion disks (e.g., see Greenhill 
et al. 1997, Moran 1998) and by the finding that in some 
contexts, especially in the dynamics of protostellar disks, 
there are empirical indications that the amount of mass 



in the disk is large (see Hillenbrand et al. 1992, Drimmel 
1996). From the theoretical point of view, to address the 
dynamics of a self-gravitating accretion disk is highly at- 
tractive especially because some global regulation aspects 
that are emerging as interesting and relevant from various 



lines of thought (see also Coppi 198C) should find a nat- 
ural manifestation in systems where the long-range forces 
are bound to generate an inherently global behavior. 

The basic framework for the construction of models 
of accretion disks is traditionally "asymmetric". While 
the momentum transport equations are generally replaced 
by a physically based prescription that bypasses our ig- 
norance on the detailed mechanisms that are involved 



(Shakura & Sunyaev 1973; see also the recent modifica- 



tion proposed by Heyvaerts et al. 1996 and Bardou et al. 
|1998 ), the energy transport equations are usually kept in 
their ideal form, either following a detailed inclusion of 
the radiation transport across the disk (e.g., see Shakura 
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& Sunyaev 1973; Bardou et al. 1998) or by invoking some 
ideal "equation of state" and by arguing that the energy 
dissipated by viscosity is partly redistributed in the disk 



(see Narayan & Yi 1994 and many following articles). In a 
situation where we lack strong empirical constraints (such 
as those that might be available in the laboratory) on 
the detailed mechanisms at the basis of the various trans- 
port processes involved, one might try to explore models 
where self-gravity, much like in galaxy disks, plays a ma- 
jor role at all the three levels mentioned above, taking 
the view that both momentum and energy transport equa- 
tions should be handled heuristically. This is inde ed th e 
step taken in an earlier exploratory paper (Bertin 1997), 
where th e visc osity prescription suggested by Shakura & 
Sunyaev ( 1973| ) is retained, and, on a similar footing, the 
energy transport equations are replaced by a physical con- 
dition of self-regulation, related to marginal Jeans stabil- 
ity, as suggested by the dynamics of galaxy disks (see Sect. 
3). ^ 

In this paper we develop this latter point of view by 
describing the entire parameter space available for self- 
regulated accretion disks (Sects. 2 and 4), and by providing 
two non-trivial extensions of the model. The first exten- 
sion is that of a "bimodal" disk, where full matching with 
an inner "standard" Keplerian accretion disk is obtained 
(see Sect. 5.1). The second extension is the construction 
of self-regulated accretion disks embedded in a spheri- 
cal diffuse "halo", which may find application to the ex- 
tended disk associated with AGN's. The analysis is further 
strengthened by a careful discussion of the vertical struc- 
ture of the disk, in such a way that the transition from 
self-gravity dominated to non-gravitating disks is covered 
uniformly (Appendix A). 

2. Self-regulated accretion disks 

2.1. Basic equations and parameter space 

We consider a steady, axisymmetric, geometrically thin 
accretion disk, with constant mass and angular momen- 
tum accretion rates, which we call M and J, respectively. 
The conservation laws for mass and angular momentum, 
in cylindrical coordinates, are: 



M = —2TTrau, 



J = Mr'^n + 2T:uar' 



(1) 
(2) 



where a is the surface density of the disk, u is the radial 
velocity, is the angular velocity, and v is the viscosity 
coefficient. Note that we take M > in the case of infiow. 
Usually on the right hand side of Eq. (||) the two terms 
have opposite signs; the first term describes convection 
of angular momentum, while the second (negative in the 
common situation where dO, /dr < 0) accounts for the an- 
gular momentum transport due to viscous torques. Thus 
J > corresponds to a net angular momentum influx. 



For a cool, slowly accreting disk, the radial balance of 
forces requires: 



ld$„ 



r dr 



(3) 



where Mi, is the mass of the central object and $o- is the 
disk contribution to the gravitational potential. The radial 
gravitational field generated by the disk can be written as: 



dr 
r 



(r, z) = - 



z 

rr' 



fc2 



ka{r')dr' , 



(4) 



where E{k) and K{k) are complete elliptic integrals of the 
first kind, and k^ = Arr' /[{r + r'Y + z^] (see Gradshteyn 
& Ryzhik |1980| ). The field d$^/d r in the equatorial plane 
is obtained by taking the limit z — s- 0. Here we have pre- 
ferred to refer directly to the field (and not to the po- 
tential, as was done by Bertin 1997) and to the formulae 
applicable to the case z ^ Q; this representation is more 
convenient in view of the numerical investigations that we 
have in mind (see Sects. 4 and 5). 

For the viscosity we follow the standard prescription 
(Shakura 



Sunyaev 1973) 



V = ach, 



(5) 



where c is the effective thermal velocity and h the half- 
thickness of the disk. Here a is a dimensionless parameter, 
with < a < 1; in the following we take it to be a con- 
stant. 

In the case of dominant disk self-gravity, one may 
adopt the requirement of hydrostatic equilibrium in the 
z direction for a self-gravitating slab model, which gives 



h = 



nGa' 



(6) 



A more refined analysis of the vertical equilibrium is pro- 
vided in Appendix A. Substituting Eq. (^ into Eq. (|^) we 
obtain ua — {a/TrG)c^ , which inserted in Eq. 



yields: 



GJ = GMr^n + 2ac^r 



(7) 



In order to close the set of equations (Eq. (||), Eq. (^, 
and Eq. (|^)), we need an additional relation. In standard 
studies this is provided by an energy transport equation 
(see, for exam ple, P ringle 1981; Narayan & Popham 1993 ; 
Narayan & Yi 1994). Here we consider an alternative sce- 
nario (Bertin 1997) where the energy equation is replaced 
by the requirement of marginal Jeans stability: 



irGa 



= Q«1, 



(8) 



where k is the epicyclic frequency. This assumes the pres- 
ence of a suitable self-regulation mechanism (see Sect. 3 
below). For a given value of a and Q, we now have a com- 
plete set of equations, depending on the three parameters 
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M^,, M, and J. We recall that the parameter cK/nGa is 
the fluid analogue of the axisymmetric stability parameter 



introduced and described by Toomre (1964) in his analysis 
of the stability of a stellar disk. 

2.2. The self-similar pure disk solution 

In the special case of J = and = 0, the problem is 
solved by the self-similar solution with flat rotation curve 



(Bertin 1997), characterized by: 
2nG<7r ^r^fl^ ^V^ ^ const. 



(9) 



The other properties of the disk are specifled by: c — 
(_GA//2a)i/3, V = (2V2/Q)c, u = -{aQ^/A)c, h/r = 
Q^/A. Note that Eg. (|^) describes the self-similar disks 
found by Mestel ( 1963[ ) in the non-accreting case. 

The same solution is also valid asymptotically at large 
radii even when M^, and J do not vanish. In fact, the effect 
of a non-zero M* should be unimportant for 

r > = 2GA/,(g/4)2(GAf/2a)-2/3 

while the effects of a non-zero J should be unimportant 
for 

r > rj = V2(| j|/M)(Q/4)(GM/2a)^i/3 

. Here the definitions of r^ and rj are slightly different 
from those adopted earlier (Bertin 1997). 



2.3. The iteration scheme for the general case 

For a given value of a and Q and for a specified accretion 
rate M, the general case is a problem with a well-defined 
physical lengthscale rg and one dimensionless parameter 
^ = sgn{ J)rj/rg. In the following we wish to explore the 
properties of the mathematical solutions in the entire pa- 
rameter space available (in particular, for both positive 
and negative values of ^). The knowledge of these solutions 
is a prerequisite for a discussion of any specific astrophys- 
ical application. At a later stage, the physical problem 
under investigation is expected to restrict the relevant pa- 
rameter space. 

In order to calculate such one-parameter family of solu- 
tions, we start by writing the relevant equations in dimen- 
sionless form, in such a way that the self-similar solution 
described above is easily recognized. We thus introduce 
three deviation functions p, (j), and x, which are related 
to the dimensional disk density a, rotation curve V , and 
effective thermal speed c by the following definitions: 



y2 

2ttG(j = -2-(1 
r 



P), 



= Vr 



2 J,2 



( Q 



\2^/2 



X, 



(10) 



(11) 



(12) 



with = (8/0^)(GA//2a)2/3. Clearly the functions 
(1 -I- p), (jy^ , and x s-re all positive definite; the self-similar 
solution corresponds top = 0, (/) = l,x=l. 

In terms of the dimensionless radial coordinate x = 
r/js, the basic set of equations (corresponding to Eqs.(||), 
(|), (0), and (1)) thus become: 



X 



dx'- 



:X 



4 1 - fc2 



\/xx 

x' x S . , 
7 + — Eik) 



X 



X' 



XX' 



p{x'), 



X(j) 



1/3 



1 - 



d In ( 
din: 



-1/3 



p{x) = x4>\ + 



din a; 



(13) 



(14) 



(15) 



(16) 



with P = 4:xx'/[{x + x'f + S'^], in the limit ^ ^ 0. 

In the case of outward angular momentum flux (i^ < 0) , 
the above equations can be readily solved by iteration in 
the following way. An initial seed solution (/jq — y/l + l/x 
is inserted in Eqs. (|l5|) and ([l6|), thus leading to a first 
approximation to the density deviation po{x\ This is in- 
serted in Eq. (Q), then producing (via Eq. (p^) a new ex- 
pression for the rotation curve deviation 4ii{x). Typically 
three or four iterations are sufficient to reach a satisfactory 
convergence. 

For the case when the net angular momentum flux is 
inwards (^ > 0), Eq. (|l^) shows that we may run into 
a difficulty at small radii, where the quantity (1 — ^/xcf)) 
changes sign. Here we may proceed by analogy with stan- 
dard studies (e.g., see Pringle 1981 ), by restricting our 



analysis to the outer disk defined by x > Xin, with the 
condition Xin4){xin) = J = -^^^in^(^"i)- This is 

taken to occur at a point where the angular velocity ^{r) 
reaches a maximum (so that dln^/dlnx — 1 and x niay 
remain finite); such maximum is identified as a location, 
in the vicinity of the surface of the central object, to which 
the accretion disk is imagined to be connected by a rela- 
tively narrow boundary layer. Of course, in the boundary 
layer the physical processes will be different from those de- 
scribed in our model (for example, there will be effects due 
to pressure gradients and the disk will be thick). Thus we 
will not be able to follow the associated (inwards) decline 
of the angular velocity away from the Keplerian profile. 
As a result, in our calculation we let the effective ther- 
mal speed vanish at rin] this unphysical behavior would 
be removed when one describes the boundary layer with 
more realistic physical conditions (see Popham & Narayan 
|1991D . 
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3. Self-regulation in self-gravitating accretion 
disks 

Before proceeding further, we should make here a digres- 
sion and try to better explain the physical justification 
at the basis of the regulation prescription set by Eq. 
which is the distinctive property of the class of models ad- 
dressed in this paper. This discussion expands and clarifies 
a brief description provided earlier (Bertin |1997 ). 



3.1. The mechanism 

We start by noting that, when self-gravity dominates, pro- 
cesses associated with the Jeans instability basically de- 
termine the relevant scales of inhomogeneity of the system 
under consideration. On the other hand, it is well known 



(Toomre 1964 ) that in the plane of a thin, self-gravitating, 
rotating disk the tendency toward collapse driven by self- 
gravity, which is naturally expected to be balanced locally 
only at small wavelengths by pressure forces, can be fully 
healed by rotation also at small wavenumbcrs. As a result, 
if the disk is sufhciently warm, it can be locally stable 
against all axisymmctric disturbances. The condition of 
marginal stability thus ensured is usually cast in the form 
(5 = 1, where Q is proportional to the effective thermal 
speed of the disk. 

Therefore, in galactic dynamics it has long been rec- 
ognized that an initially cold disk is subject to rapid evo- 
lution, on the dynamical timescale. For a stellar disk, the 
Jeans instability induces fast overall heating of the disk 
up to levels of the effective thermal speed for which the 
instability is removed, basically following the above crite- 
rion (Q « 1). The heating rate is thus very sensitive to the 
instantaneous value of Q. Therefore, even in the absence 
of dissipation, for relatively low values of Q the collective 
Jeans instability provides a mechanism to stir the system 
and to heat it. However, a coUisionless non-dissipative disk 
would be unable to evolve in the opposite direction, i.e. to 
cool if initially hot {Q > 1), and, because of other residual 
processes (such as tidal interactions), may still be subject 
to some perennial heating even when Jeans instability is 
ineffective. In turn, if efficient dissipation is available (such 
as that associated with the interstellar medium in galaxy 
disks), a competing mechanism can cool an initially hot 
disk down toward conditions of marginal Jeans instability. 

These important dynamical ingredients have been 
studied from various points of view, giving rise to inter- 
esting scenarios, also by means of numerical experiments 



(see Miller et al. |1970[ Quir k |1971| , Q uirk fc Tinsley |1973 
Sellwoood & Carlberg 1984, Ostriker 1985), so that even- 



tually the possibility of an actual self-regulation in self- 
gravitating disks has been formulated and explored. Ac- 
cordingly, the two competitive mechanisms, of dynamical 
heating and dissipative cooling, can set up a kind of dy- 
namical thermostat. In particular, it has been shown that 
self-regulation is very important in determining the condi- 



tions for the establishment of spiral structure in galaxies 
(Bertin and Lin 1996; and references therein). The gen- 
eral concept has also found successful application to the 
dynamics of the interstellar medium, even when dynam- 
ically decoupled from the stellar component, in relation 
to the interpretation of the observed star formation rates 



(Quirk 1972; Kennicutt 1989). Note that the most delicate 
aspect of the regulation process is the cooling mechanism; 
in galaxy disks, an important contribution to cooling is 
thus provided by cloud-cloud inelastic collisions that dis- 
sipate energy on a fast timescale. 

A final remark is in order. As shown by the case of 
galaxy disks, the dynamics involved may be extremely 
complex, so that there is little hope to provide a simple 
"ideal" energy equation for the description of a coupled 
disk of stars (of different populations) and gas (in differ- 
ent forms and phases). This point has an additional impor- 
tant consequence. If we try to describe such an inherently 
complex system by means of an idealized one-component 
model, we should be ready to introduce the use of some 
effective quantities (in particular, of an effective thermal 
speed; see also the use of the term "effective" after Eq. ^ 
and at the beginning of this subsection). 

3.2. Viability of the mechanism for some accretion disks 

It should be emphasized that the main strength of the 
self-regulation scenario described above is rooted in semi- 
empirical arguments. The actual data from many galaxy 
disks (among which the Milky Way Galaxy) and from 
planetary rings show that conditions of marginal stabil- 
ity (Q « 1) often occur and can indeed be established in 
systems subject to complex dynamical processes. For ac- 
cretion disks, there may already be some empirical indica- 
tions in this direction, to the extent that the application 
of standard models to some observed systems points to 
very cold disks, with Q well below unity (e.g., see Kumar 
1999, as briefly mentioned in Sect. 6). Additional clues 
in the same direction also derive from numerical experi- 
ments; in the context of protostellar disks, numerical sim- 
ulations show that disks formed under conditions where 
self-gravity is important include wide regi ons c haracter- 
ized by a constant Q-profile (Pickett et al. 1997 ). 

An accretion disk may be subject to efficient cooling by 
a variety of mechanisms, depending on the physical con- 
ditions that characterize the specific astrophysical system 
under consideration. From this point of view, the cooling 
necessary for the establishment of self-regulation may oc- 
cur efficiently already via the radiative processes included 
in "standard" models (see the general discussion of the 
outer disk by Bardou et al. 



1998, which we will summa- 



rize in Sect. 5.1). In reality, the dynamics of matter slowly 
accreting in a disk can be significantly more complex. Cold 
systems, such as a protogalactic disk, a protostellar disk, 
or the outer parts of the disk in an AGN, may have a com- 
posite and complex structure. They may include dust, gas 
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clouds, and other particulate objects with a whole variety 
of sizes and "temperatures" . Much like for the HI compo- 
nent of the interstellar medium, the main contribution to 
the effective temperature of the disk might be from the 
turbulent speed of an otherwise cold medium. On the one 
hand, for these systems it may be hard or even impossi- 
ble to write out a simple "ideal" energy transport equa- 
tion. On the other hand, such a complex environment is 
likely to possess all the desired cooling and heating mech- 
anisms that cooperate in self-regulation. In this respect, 
one is thus encouraged to bypass the problem of defin- 
ing a representative set of equations for energy transport, 
and to use instead the semiempirical prescription of Eq. 

Somewhat in a similar way, our inability to derive 
from first principles a satisfactory set of equations for mo- 
mentum transport is often taken to justify the adoption 
of the a-prescription of Eq. (||) . These phenomenological 
prescriptions have several limitations, but may still work 
as a useful guide to our efforts and provide interesting 
models to be compared with the observations. 

To be sure, some types of accretion disk, or some re- 
gions inside accretion disks (for example, very close to the 
center; see Sect. 5.1), may lack the physical ingredients 
invoked above. In fact, there is no reason to claim that 
self-gravity must always be important. Therefore, we will 
study the structure of self-regulated accretion disks, as a 
viable class of astrophysical systems, while we do recog- 
nize that warmer, non-regulated disks may exist and are 
likely to be basically free from the effects associated with 
the self-gravity of the disk. 

3. 3. The impact of Q on momentum and energy transport 
The self-regulation mechanism has been demonstrated by 



ization of the actual system that we are addressing. Given 
the indications of several dynamical studies (in addition 
to those of Lin & Pringle, see, for example, Laughlin & 
Bodenheimer 1994, Laughlin & Rozyczka 1996), a more 



considering a simplified set of equations (Bertin 1991) 
where efficient cooling is included and the role of self- 
gravity is modeled by means of a heating term with an an- 
alytic expression (inversely proportional to a high power of 
Q) meant to incorporate the results of dynamical studies 
that show that heating is indeed very sensitive to the value 
of Q. The main features of this formula, with its threshold 
at Q w 1, aimed at representing the "thermal evolution" 
of the disk, are somewhat analogous to the heuristic char- 
acterization of the v iscosity dependence on Q adopted by 
Lin & Pringle ( 1990 ) in the parallel problem of construct- 
ing the momentum transport equations when self-gravity 
is important. 

In our discussion of self-regulated accretion disks, we 
actually have no doubt that self-gravity is likely to have 
an important impact on viscosity, and this is still tacitly 
incorporated in the a prescription. This impact is even 
more obvious if one recalls that a self-gravitating disk 
can be subject to non-axisymmetric instabilities, which 
are bound to contribute significantly to angular momen- 
tum transport. Our class of axisymmetric, steady-state 
accretion models represents only one approximate ideal- 



complete analysis should thus include one further relation 
between a and Q. In general, this might practically re- 
quire that the phenomenological prescription (^ be used 
with a parameter a varying with radius. In reality, we be- 
lieve that the proper way to include the relevant physical 
effects, especially those associated with non-axisymmetric 
instabilities, would be through some global constraint. Un- 
til such global description remains not available, the as- 
sumption of a free, constant a may provide a first approx- 
imation, best applicable when Q is self-regulated. This 
choice can be physically consistent a posteriori, at least 
for the self-similar solution of Sect. 2.2. 



4. Properties of models with a central point mass 

After this digression, we can now proceed with the analysis 
of the problem as formulated in Sect. 2. In the presence of 
a central point mass (M^ ^ 0) we expect the lengthscale 
rg to mark the transition from a Keplerian disk to a fully 
self-gravitating disk with flat rotation curve. Surprisingly, 
in our class of self-regulated accretion disks the role of 
the disk self-gravity turns out to be significant all the way 
down to the center. 



J^.l. Rotation curves 

In Fig. we illustrate the behavior of the rotation curves 
in our class of models for several values of the angular 
momentum flux parameter For comparison, we also 
show the Keplerian curve Vk that is obtained by setting 
1-1-/9 = 0. We note that the difference from the Keplerian 
decline is significant even at r = 0{rs)- For example, for 
the ^ = case we find {V - Vk)/Vk ~ 100% at r = 2r^. 



4-2. Effective thermal speed 

Disks where the angular momentum is transported out- 
wards (^ < 0) tend to develop a warmer core, while the 
opposite trend occurs for disks where the angular momen- 
tum is carried inwards {£_ > 0). This is shown in Fig. || 
(Fig. H and the following Fig. |3[ Fig. ^, and Fig. || are 
shown in logarithmic scale to better bring out the behav- 
ior in the inner parts of the disk). Note that in the case 
of inward angular momentum flux (^ > 0) the effective 
thermal speed need not vanish at the inner edge of the 
disk (as it does in our models) where a boundary layer 
is expected to be generated (see discussion at the end of 
Section 2.3). 
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Fig. 1. Top panel: Rotation curve of the disk in the presence 
of a central point mass for different values of the angular mo- 
mentum flux parameter. Bottom panel: Same as top panel, 
but in logarithmic scale. 

4-.3. Role of the disk self-gravity close to the central point 
mass 

There are several quantities directly related to the disk 
density distribution in the disk that allow us to character- 
ize the role of the disk self-gravity. 

The most natural quantity to consider is the ratio of 
the mass of the disk to that of the central object. Ob- 
viously, for our non-truncated models this quantity is 




lri(r/rj 

Fig. 2. Equivalent thermal speed of the disk. The two cases 
> and ^ < 0) show opposite behavior in the inner disk. 
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Fig. 3. Cumulative mass of the disk relative to that of the 
central object. 

meaningful only when referred to a given radius. Fig. |^ 
shows how rapidly in radius the system becomes domi- 
nated by the mass of the disk. Note that, in any case, 
Md^sfcW/M^^Oforr^O. 

In galactic dynamics the local disk self-gravity is usu- 
ally measured in terms of the parameter e = irGa/rK^. 
The fully self-gravitating self-similar disk (with flat rota- 
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Fig. 4. Local self-gravity of the disk, as measured by e = 
-kGo/th^, for different values of the angular momentum flux 
parameter. 
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Fig. 5. Ratio of the vertical scaleheight ft* = cr/Vx to the 
thickness of the disk for various models. 

we have also considered models based on the improved 
prescription for the disk thickness (see Appendix A): 



tion curve) is characterized by e = 1/4. The profile of this 
parameter (see Fig. ^) confirms that indeed, close to the 
center, the influence of the central mass becomes stronger 
and stronger. 

Given the behavior of the profiles Mdiskir) 1^* and 
e(r), one might conclude that the innermost disk should 
be treated as a standard Keplerian accretion disk. This 
conclusion is contradicted by the following argument. In 
setting up the equations of our models, we have taken 
the vertical equilibrium to be dominated by the disk self- 
gravity (see Eq. (^). If the innermost disk were fully Ke- 
plerian, at small radii the vertical scaleheight = cr IVk 
should become much smaller than the thickness h asso- 
ciated with our models. Instead, a plot of the ratio 
only shows that the two scales become comparable to each 
other (see Fig. ||), thus demonstrating that the influence 
of the disk self-gravity is significant all the way to the 
center. This, of course, reflects our choice of imposing the 
self-regulation prescription (Eq. (||)) at all radii; in the 
next section we will show that the disk can indeed make 
a true transition to a Keplerian disk if the self-regulation 
prescription is suitably relaxed. 

In view of the above considerations, in order to check 
that no major consequences arise from the use of a vertical 
equilibrium prescription only partly justified for r < Ts, 



TTG<jAQ^2nyK^-i) ^^^^ 
I 8 , f2n^ V 

Note that in the limit 2^ j — > 1, applicable to the self- 
similar disk, this prescription reduces to Eq. (|^). In the 
innermost region — > Vk/t, — > 1, so that h^^/h 
(4gA)/[v/l + 8Q77r-l]. 

In Fig. H we show the rotation curve of the improved 
model, compared to that of the original one, for the £, = 
case. This plot (along with similar results for the other rel- 
evant physical quantities of the disk) shows that, qualita- 
tively, the more refined vertical analysis leaves the models 
basically unchanged. 

5. Extensions 

5.1. Matching with an inner Keplerian, 
non- self -regulated accretion disk 

If we take an astrophysical situation (such as an AGN or 
a protostar) with specific physical conditions, it is likely 
that the arguments that support the adoption of the self- 
regulation constraint, followed in this paper, fail outside a 
well-defined radial range, either at small or at large radii. 
For example, we may recall that in the context of the dy- 
namics of spiral galaxies the relevant Q profile is argued 
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Fig. 6. Rotation curves for two different treatments of the 
vertical structure for the case ^ = 0, Q = l.The solid line 
represents the improved model, based on Eq. (|l7|); the dotted 
line is the original model, based on Eq. 



to be flat in the outer disk, but is thought to increase in- 
wards inside a circle of radius rg often identified as the 
radial scale of influence of the bulge; even in the absence 
of a bulge, the central parts of the disk a re tho ught to 
be generally hotter (e.g., see Bertin & Lin, 199£). In our 
context, we may then imagine an accretion disk where at 
small radii self-regulation fails, Q becomes large, and, cor- 
respondingly, the disk self-gravity ceases to be important. 

This discussion suggests that it should be interesting 
to explore the possibility where Eq. (|^) is replaced by a 
condition of the form cK/nGa — Q{r), with a Q profile 
that decreases monotonically with radius and reaches the 
self-regulation value Q only beyond a certain radial scale 
rg. What would be the impact of such a choice on the 
structure of the accretion disk? What might be the physi- 
cal arguments leading to the justification of such a profile? 
In particular, what would set the scale rg for given values 
of the parameters a, M^,, M, and J? Note again that we 
are reversing the standard point of view, whereby one may 
ask what is the Q profile for an accretion disk, based on 
the structure calculated from a given choice of the energy 
equations. 

Imposing, as we are going to do, a given profile Q(r) 
may, at first sight, appear to be arbitrary. In reality, the 
freedom in the choice of the profile allows us to test quanti- 
tatively how the dynamical characteristics V(r), Mdiskir), 
e(r), and h{r) of the disk change when the self-regulation 
constraint is partially relaxed, in the inner disk, in a vari- 
ety of ways. It is up to us to test the different possibilities 



(which may correspond to completely different sets of en- 
ergy balance equations in the inner disk) that might be 
considered. For our purposes, since we wish to study the 
deviations from the standard "Keplerian" case, we only 
need to take into account that the relevant physical pro- 
cesses match so that in the outer disk self-regulation is 
enforced. Note that in the transition region where match- 
ing between the inner and the outer disk occurs, for rea- 
sons expressed in Sect. 3, it may be practically impossible 
to define, from first principles, a satisfactory set of en- 
ergy equations able to include all the desired radiation 
processes and the non-linear effects associated with Jeans 
instability. 

This procedure draws considerable support from a re- 
cent analysis of "standard" disks (Bardou et al. 1998 ) 
aimed at detecting evidence for the importance of disk 
self-gravity in the outer disk. Based on an extension of 
the "standard" a-disks (characterized by Kramers' opac- 
ity and by neglect of radiation pressure) , the effects of the 
disk self-gravity have been here incorporated by means 
of an improved thickness prescription (somewhat in the 
spirit of our Appendix A) and of a modified viscosity pre- 
scription, but the (Keplerian) profile is left unaltered. 
Therefore, this study is ideally suited to describe the con- 
ditions of our inner disk, as we intend to partially relax 
the self-regulation requirement. A very important result 
of the analysis by Bardou et al. ( |1998| ) is that their "stan- 
dard" description breaks down beyond a radius rg, well 
inside which the local stability parameter behaves approx- 
imately as Q ^ J.-9/8. might have been anticipated, 
the location where the standard model breaks down coin- 
cides with the location where Q becomes approximately 
eq ual to unity. In conclusion, the analysis by Bardou et. 
al ( 1998 ) encourages us to consider the following choice of 
Q profile 



TTGa 



= Q{l + ir/rQr'/'eM~ir/rQ)]} 



(18) 



to be used instead of Eq. (||). (The formula is meant to 
be used as a semi-empirical tool; one should keep in mind 
that the exact form of the Q profile will be determined 
by the detailed energy processes occurring in the inner 
disk and by the progressively important role of collective 
instabilities.) If we express the scale rg found in that study 
in terms of our scale r^, we find 



^ -0.8-10-4- 



-2/45 / 



g2 V io^Mq 



-2/3 



M 



8/45 



lOMQ/yr 



.(19) 



Note that the ratio rg /rs decreases while M^, increases, 
and that its value is only weakly dependent on a and on 
A'l. For parameters typical of an AGN, self-regulation may 
thus be ensured very far in, while an extrapolation of the 
above recipe to parameters typical of a disk surrounding 
a T Tauri star suggests that, for these latter objects, rg 
should become 0{rs). 
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Fig. 7. Rotation curves of partially self-regulated models with 
the Q-profile given in Eq. (1^), for different values of rg. The 
thickness prescription used is that of Eq. (0). 



Fig. 8. Rotation curves of models with a diffuse halo, with 
5 = 0, and core radius equal to (a;o = 1), for different values 
of f. 



Some examples of partially self-regulated models com- 
puted on the basis of Eq. (18) (see Fig. ^ show how an 
outer disk dominated by self-gravity can match in detail 
with an inner standard Keplerian accretion disk. 

5. 2. The effect of a diffuse "halo " 

So far we have considered the case where the mass is all 
distributed in a disk (either at the center, as a point mass, 
or in diffuse form) . In view of possible applications to AGN 
configurations or to the general galactic context, it is im- 
portant to consider an extension of the models to the case 
where part of the gravitational field is determined by a 
diffuse spherical component, which we will call halo (even 
if it may just correspond to the central region of an ellip- 
tical galaxy). This will lead to rotation curves otherwise 
not accessible by our models. On the other hand, it is eas- 
ily recognized that this natural extension is going to leave 
the slow density decline of the disk unaltered. Therefore, 
if we are interested in producing models with finite total 
mass, we should be ready to impose an outer truncation 
radius or, which might effectively be equivalent, to relax 
the self-regulation prescription in the outermost disk. 

We have thus considered a set of models where the field 
external to the disk is produced by the joint contribution 
of a central point mass and of a halo (which, for simplicity, 
we take to be spherical). In view of the case of a disk 
embedded in an elliptical galaxy, we have modeled the halo 
as approximately isothermal, with a finite core radius. In 



this case the dimensionless equation giving the rotation 
curve (Eq. ([l^)) is modified as follows: 



(20) 



We see that now the equations depend on two additional 
parameters: giving the relative strength of the external 
field, and xq, which measures the size of the core radius. 
In this case it is easy to demonstrate that at large radii 
the density deviation p approaches if / <C 1, and / if 
/»!. 

In Fig. ^ we show examples of the rotation curve of 
models with a diffuse halo, for the case ^ = 0, xq = 1. For 
the vertical structure, we have referred to the improved 
vertical prescription of Eq. (A. 9), with fl^^t oc f'^/{xQ + 



'-)■ 



5.3. Disks with an outer truncation radius 

Self-regulated accretion disks with finite mass can be eas- 
ily constructed by imposing the existence of an outer trun- 
cation radius. Either the study of the collapse of a gas 
cloud with finite mass or the consideration of the physical 
conditions in the outer parts of some astrophysical objects 
will naturally bring us to address such models. 

6. Discussion 

In this paper we have developed a framework for the con- 
struction of a class of self- gravitating accretion disks, in- 
dependently of the specific conditions characterizing the 
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astrophysical systems where the accretion disk paradigm 
apphes. In principle, we might even speculate that the 
above framework could be the starting point to describe 
some stages of the formation of protogalactic disks (e.g. 
see Fall & Efstathiou 198C) or the establishment of some 
regular extended HI disks in eUiptical galaxies (see Mor- 
ganti et al. 1997). In practice, our class of models appears 



to be best applicable to two categories of astrophysical 
objects, i.e. to Active Galactic Nuclei and to protostellar 
disks. Concrete applications, which will require a detailed 
consideration of the available observational constraints, 
are beyond the scope of the present paper. Still we would 
like to make a few comments that should illustrate why the 
models appear to be promising for the above categories of 
objects. 

Recently, from radio maser emission, it has been possi- 
ble to obtain accurate measurements of the rotation of the 
disk in the central parts of a few Active Galactic Nuclei; 
these include NGC 4258, NGC 1068, Circinus, NGC 4945, 



NGC 3079, and NGC 1386 (see Moran, |l99§). These mea- 
surements have shown that, although in some cases (as for 
NGC 4258, see Miyoshi et al. 1995) the rotation curve is 



Keplerian to a high degree of approximation, there may be 
significant deviations from the r~^l'^ profile. For example, 
in NGC 1068 the rotation velocity at distances « 1 pc from 
the "center" turns out to decline as r~° '^^ (see Greenhill 



et al. 1997; Kumar 1999). The general conclusion is that 
rather extended disks, which may contain large amounts 
of mass, are present; indeed one possible reason for the 
deviation from the Keplerian decline has been identified 
in the influence of the gravitational field contributed by 
the accretion disk itself. At the same time, the applica- 
tion of "standard" accretion disk models to some of these 
objects has led to finding values of the Q parameter well 



below unity (for NGC 1068, see Kumar |1999D . AU these 
are clues that show that a model where the role of the disk 
self-gravity is fully incorporated is called for. Note that if 
we adopt the numbers suggested by the data for NGC 
1068 we would find Ts « 25 pc and tq k, 1.6 lO^'^r^; then, 
it is curious to find that a self-regulated disk with these 
characteristics would have, for ^ = —5, a gradient of the 
rotation curve at r w 1 pc compatible with V ^ 7.-0.37^ 

On a completely different mass scale, if we refer to 
the case of protostellar disks, typically quoted parameters 
are M ? ^ lO^^Mo/j/r and Ah « O.5M0 (Hartmann et 
al. 1998 ). Under these circumstances we find Vg « 1000 
AU and rq « 0.8rs. Interestingly, protostellar disks have 
been observed to extend out to a radius from sa 100 AU to 



1000 AU (Dutrey et al. |1996| ; McCaughrean & O'Dell 

1996| ). As for the case of AGN's, a check on the values 



of the temperatures anticipated on the basis of the effec- 
tive thermal speeds predicted by the self-regulated models 
shows that the numbers fall reasonably within the range 
suggested by the observations. 

Finally, it is interesting to note that studies of the 
spherical, inviscid collapse of a molecular cloud, imag- 



ined to eventually generate a system composed of a proto- 
star and a circumstellar disk, lead to mass accretion rates 



M (X (Hunter 1977, Shu 1977), curiously analogous to 
those of our completely different accretion scenario. This 
coincidence may offer interesting clues for modeling the 
entire process where a cloud, starting out in an initially 
spherical collapse, ends up in a protostellar accretion disk. 
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Appendix A: Modeling the vertical density profile 

The vertical structure of the disk is generally derived by 
imposing hydrostatic equilibrium in the z-direction un- 
der the assumption of isothermality. The relevant Pois- 
son equation and hydrostatic equilibrium condition for an 
axisymmetric configuration can be written as (here the 
symbol p denotes volume mass density and should not be 
confused with the density deviation function of the main 
text): 



-TT r— + — y = 47rGp + AirGpexu 
r or \ or J oz^ 



(? dp 
p dz 



dz ' 



(A.l) 



(A.2) 



where we have taken polar cylindrical coordinates. For 
simplicity, we assume that the external field associated 
with pext is spherical. [In the class of models studied in 
Sect. H and Sect. ^ we have taken the case where pext = 
Mj,(5(x).] From the definitions of Vl and of k, it is re adily 
shown that, close to the equatorial plane, Eq. ( A.l ) can 
be rewritten as: 



92$ 

dz"^ 



= At:Gp + AnGpext + - k'^ 



AttGp + (2f72 - k2 



(A.3) 
(A.4) 



The latter expression follows from the exact relation 
AirGpext — {k^ — ^'^)ext applicable to a spherical density 
distribution. The expression (217^ _ k^)^. indicates that 
the quantity in parentheses should be calculated based 
on the radial field given by Eq. (Q). For a relatively thin 
disk the quantity (2il2 — k^)^. -f fl'^^^ may be taken to 
be nearly independent of z, s o th at the gravitational field 
obtained by integrating Eq. (A.4) includes a term that is 
approximately linear in z. The result can be inserted in 
the right hand side of Eq. (A.2), which becomes: 



dV. 

dz2 



dz 



27rG 



(2f^2 - k2)^ 



02 



4,(A.5) 



where we have introduced the integrated density cr^ — 
2 p{z')dz' . Note that the term ^^^tZ is just the vertical 
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component of the spherical external field, as might have 
been anticipated. If po = p{z = 0) is the density on the 
equatorial plane, we can define two scales, S = yJIpaj-KGc 
and a = c/y/2TrGpo, and the natural dimensionless pa- 
rameter A = [(2r^2 - K'^)a + f2ext]/47rG'/9o, so that Eq. 



(A. 5) becomes: 



dC ^ 



(A.6) 



to be solved under the conditions y{0) — 0, y'{0) — 1; 
here we have used the dimensionless variables y — (Tz/S 
and C = z/H. For each value of A one can compute the 
desired density profile and the associated surface density 
a = az{z = oo). One can then introduce a scaleheight h 
such that a = 2pQh, and the value of h can be computed 
directly from the relation h — Hy{z = oo). 

The standard non self-gravitating model, where the 
contribution of the gas density p to the vertical gravita- 
tional field is negligible, corresponds to the limit equation 
y" — —2A(y'. Usually the analysis is carried out with 
the additional approximation A « rtl^^/AwGpQ, but this 
is not needed; by retaining the contribution (2fl^ — k}^u 
one can thus keep track of the effects associated with the 
rest of the disk mass distribution, which may be signif- 
icant even where the local disk density po is small. In 
any case, such limit equation yields the well known Gaus- 
sian profile p = Po exp(— with /ij — H'^/2A 



the standard expression hi, « c/^k used in the so-called 
Keplerian limit. Note that hi, is slightly different from the 
scaleheight h defined above, since h = yjir /2hi,. 

The limit of the homogeneous fully self-gravitating 
slab corresponds to the equation y" = —2yy', so that 



the vertical density profile is given by (Spitzer 1942) 
p — Po sech^{z/h), with h — c^/irGa; here the scale h 
is the same as that defined by the relation cr — 2poh. Note 
that even when no external (spherical) field is present , the 
solution for an axisymmetric disk obtained from Eq.( A.6 ) 
is not exactly the one derived in the homogeneous, self- 
gravitating slab, unless the rotation curve is flat. 

The density profile associated with Eq.( [A.6|) is nei- 
ther Gaussian nor sech^. For practical purposes it may 
be convenient to use a simple interpolation formula for 
the vertical scale, which is justified by the following 
description "biased" towards the Gaussian limit. If we 
start from Eq.( |A.2| ), naively expand the vertical field as 
{d<S>/dz){z) - {d^<^>/dz^),^oz, and then use Eq.(|01) to 
estimate the factor {d'^^/dz'^)z=a, we find an equation 
leading to an unrealistic Gaussian density profile with 
scaleheight h given by: 




A 



Fig. A.l. Comparis on l aetween the exact (num erica lly com- 
puted based on Eq. (A.6)) and approximate (Eq. A. 8) expres- 
sion for the disk thickness. 



The quantity po contains a dependence on h (for given 
a). This suggests the use of the following interpolation 
formula: 



here improved with respect to _ 



1 



nCa J l + AA/ 



TT 



(A.; 



with f{Al^[y{z = 
guar- 



Note that the exact calculation, from h ~ Hy{z = oo) 
gives a relation h = [c^ /TTGa)f{ 
oo)Y . The choice of the factor (4/7r) in Eq. (|A 
antees the proper limits for A — Q and for ^ — s- oo. 
The accuracy of the interpolation formula is illustrated 
in Fig. (A.l). In turn, since A depends on h and on a via 
the quantity pq, for the purposes of the present paper it 

as: 



may be convenient to reexpress Eq. (A 



■nGa 4a 



(A.9) 



with a = Q'^[{2^'^-K^)a + nl^t]/K'^. This leads to Eq. ^ 
of the main text, applicable when the external field is pro- 
duced by a central point mass. This interpolation formula 



improves on earlier analyses (Sakimoto & Coroniti 1981 



Bardou et al. 1998 ), in several respects, allowing us to 
better describe the transition between a Keplerian and a 
self-gravity dominated disk and to extend the treatment 
to the case where the external field is not just that of a 
simple point mass at the center. 



TT 

2h? 



^[1 + ^]. 
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